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We study the fate of double-Weyl fermions in three-dimensional systems in the presence of long- 
range Coulomb interactions. By employing the momentum-shell renormalization group approach, 
we find that the fixed point of noninteracting double-Weyl fermions is unstable against Coulomb 
interactions and flows to a nontrivial stable fixed point with anisotropic screening effect. Moreover, 
experimentally measurable quantities such as specific heat obtain unusual logarithmic corrections. 

We also discuss implications of our results to three-dimensional materials HgCr 2 Se 4 and SrSi 2 , 
candidate materials of hosting double-Weyl fermions. 


Introduction: While characterizing the nature of 
possible exotic quantum phases in correlated electronic 
systems with large Fermi surfaces remains a central and 
challenging issue, electronic systems with only discrete 
Fermi points have attracted increasing attention in the 
past few years[T]. With only Fermi points, new phenom¬ 
ena beyond conventional Fermi liquids are often expected 
when appropriate short-range interactions or long-range 
Coulomb interactions are taken into account. For in¬ 
stance, for systems with simple linear band crossing such 
as single-layer graphenes [2] and surface states of topo¬ 
logical insulators 013] as well as three-dimensional Weyl 
semimetals, Lorentz invariance[T| and, under certain con¬ 
ditions, space-time supersymmetryjSHSj could emerge in 
low energy and long distance. For systems with nonlinear 
dispersions around band-touching points, the interaction 
effect is often more manifest due to enhanced density of 
states. The simplest such nonlinear dispersion is the case 
of quadratic band touching (QBT)[S|. It was shown that 
for two-dimensional QBT systems spontaneous symme¬ 
try breaking already occurs at infinitesimal short-range 
interactions [SHE] . 

For QBT points in three-dimensional (3D) cubic sys¬ 
tems such as Pr 2 lr 207 , weak short-range interactions 
are simply irrelevant mainly because of vanishing den¬ 
sity of states. However, long-range components of 
Coulomb interactions are relevant at the noninteracting 
fixed point and may induce a non-Fermi liquid [Ml I15j 
or a Mott insulator [161 Hz], which is in sharp contrast 
to marginally irrelevant Coulomb interactions [TSH13] in 
Weyl semimetals . Another nontrivial band cross¬ 

ing whose dispersion is not simply linear is the double- 
Weyl point, around which spectrum disperses linearly 
only along one momentum direction but quadratically 
along the two remaining directions [50115^ . As suggested 
by its name, a double-Weyl point acts like a magnetic 
monopole in crystalline momentum space with twice the 
monopole charge of usual Weyl points. Double-Weyl 
fermions may be realized in materials breaking time- 
reversal or inversion symmetry and were predicted to 
occur in the candidate material HgCr2Se4[3n| with two 


separated double-Weyl points which are protected by the 
material’s crystalline C 4 symmetry. Weak short-range in¬ 
teractions are certainly irrelevant at the noninteracting 
fixed point of double-Weyl fermions. Nonetheless, sim¬ 
ilar to the QBT in 3D, the vanishing density of states 
for double-Weyl fermions may not screen the long-range 
Coulomb interactions effectively and one naturally asks 
if the long-range Coulomb interactions can induce exotic 
new physics in double-Weyl semimetals. 

In this paper, we investigate the fate of double-Weyl 
fermions in the presence of long-range Coulomb interac¬ 
tions, particularly motivated by the candidate double- 
Weyl semimetal material HgCr2Se4[3D] with two double- 
Weyl points. (Note that two is the minimum number of 
double-Weyl points because of the no-go theorem [33].) 
Here are the main findings. At the noninteracting fixed 
point of double-Weyl fermions, the isotropic part of the 
Coulomb interaction is marginally irrelevant while its 
anisotropy is relevant; the noninteracting fixed point is 
then not stable and the system flows into a stable non¬ 
trivial fixed point with anisotropic Coulomb screening. 
The anisotropic Coulomb screening is mainly due to the 
anisotropic dispersions of double-Weyl fermions. Com¬ 
pared with the linear-dispersing direction, the quadratic- 
dispersing directions have larger density of states which 
results in stronger screening. Specifically, the Coulomb 
interactions along the quadratic directions are screened 
to a faster-decaying ^ behavior, while it is still ^ along 
the linear-dispersing direction. At this fixed point with 
anisotropic Coulomb screening, the Coulomb interaction 
is again marginally irrelevant, which is in sharp contrast 
to the anisotropic Coulomb screening in anisotropic Weyl 
fermions at whose fixed point the Coulomb interaction is 
eventually irrelevant |34j . Because of the marginally irrel¬ 
evant Coulomb interactions at the stable fixed point, var¬ 
ious physics quantities receive logarithmical corrections 
and we explicitly analyzed the temperature dependence 
of specific heat and found exotic logarithmical corrections 
at low temperature. 

Model: We start with a simple microscopic model 
of spinless fermions (namely, spin-polarized electrons 
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close to the Fermi level of the ferromagnetic half metal 
HgCr 2 Se 4 ), which hosts two double-Weyl fermions. We 
consider the following effective model featuring two 
double-Weyl fermions: 

Ho = E ct| 2 ti(cos ki — cos k2)(Ji + t2 sin ki sin fc2(T2 

k 

+ 2 t 3 (cOsfci -|- C0sfc2 -b cosfca — m)cr3|c^, ( 1 ) 

where ct creates a spin-polarized electron with orbital 

k<7 

index a — 1,2 and tj are hopping parameters. We assume 
that the system respects C 4 rotational symmetry along 
the 2 ; axis as well as the mirror symmetry changing 
z —^ —z. The phase diagram as a function of m is shown 
in Fig. 1. We shall focus on the region 1 < to < 3 where 
two double-Weyl points are realized at Q± = {0,0, ±k’^) 
with = arccos(2 — to). Specihcally, the C 4 symmetry 
protects the double-Weyl points, while the Rz symme¬ 
try requires that the two double-Weyl points have equal 
energy. 

Because of vanishing density of states at double- 
Weyl points, long-range Coulomb interactions cannot 
be screened to short-range ones. Here, for simplicity, 
we consider instantaneous 1/r Coulomb interactions be¬ 
tween electrons, which is a good approximation to the 
retarded one since the speed of light in materials is much 
faster than the Fermi velocity of electrons. By Hubbard- 
Stratonovich transformation, we introduce a bosonic field 
to decouple four-fermion Coulomb interactions so that we 
consider the following effective action to study the fate 
of double-Weyl fermions in the presence of Coulomb in¬ 
teractions: 


S = J 

— iv{—id^)r3 


{do + ig(j))To - 


iP ■ 


1 ■ 
2 . 


{d,ci>f+{d2^f+v{d3(i,r ,( 2 ) 


whCTe Fq = Ti,ri = CTir2,r2 = Cr2T2,r3 = (T3T2, 

di{k) = ki — k 2 ,d 2 {k) = 2kik2 , Tj, and tJi are Pauli 
matrices representing valley space and orbital space, and 
tjj^ = {iIjI_ , ) where i/'i represents the low-energy elec¬ 

tron near two band-touching points Q±. Note that we 
have implicitly assumed above that the coupling be¬ 
tween different double-Weyl fermions is irrelevant when 
the separation of double-Weyl points is large compared 
to the momentum scale that we are interested in (see 
Appendix). The Coulomb interaction is proportional to 
/v. In the above action, we also introduced 77 to rep¬ 
resent the anisotropy of Coulomb interactions which is 
allowed by the crystalline symmetry of the system. As 
shown below, 77 flows to zero in the infrared hxed point 
which renders qualitatively different Coulomb interac¬ 
tions in the xy plane and along the z axis. 

Screening from RPA analysis: It is well known 
that long-range Coulomb interactions between electrons 
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FIG. 1. The schematic phase diagram of the microscopic lat¬ 
tice model in Eq. 0 as a function of the parameter m. Here, 
“3D Chern insulator” refers to a 3D insulator with finite Hall 
conductivity Quadratic band touching is realized at the 
transitions between the double-Weyl semimetal and the trial 
(or Chern) insulators. 


are screened by low-energy electrons and become short- 
ranged in materials with large Fermi surface. Never¬ 
theless, in a double-Weyl semimetal, the Fermi surface 
shrinks into points and the density of states vanishes at 
zero temperature so that the screening effect needs to 
be reconsidered. In this section, we implement a ran¬ 
dom phase approximation (RPA) calculation to inves¬ 
tigate the screening effect in a double-Weyl semimetal. 
The low-energy effective action is given by Eq. 0. In 
standard RPA calculations, fermions are integrated out 
to get electron-hole polarization, the quantum correction 
to Coulomb interactions. A straightforward calculation 
shows (see Appendix for details) that electron-hole po¬ 
larization is anisotropic, 

-IV{p) <xp\ + \pz\, (3) 

where for simplicity we have implicitly neglected insignif¬ 
icant coefficients as well as a potentially-unimportant 
logarithmic factor before p\. Since the linear part of 
electron-hole polarization is more important than the 
quadratic part in the infrared, the higher order of Pz can 
be omitted at the long-wavelength limit so that we ob¬ 
tain the screened Coulomb potential Vc {p) p 2 _|_p 2 _|_|p | ■ 

This means that the bare isotropic Coulomb interaction 
becomes highly anisotropic due to the screening by the 
double-Weyl fermions. In real space, the Coulomb po¬ 
tential along the z axis is qualitatively unchanged, i.e., 
Vc{0,z) oc 1 ^, while in xy plane, the Coulomb potential 
behaves as Vc(^-LiO) oc namely, it decays qualita- 
tively faster due to the screening effects of the low-lying 
quadratic fermions in the xy plane. 

Renormalization group analysis: In this section, 
we implement a momentum-shell renormalization group 
(RG) analysis to show that rather than being screened 
all together, Coulomb interaction is marginally irrelevant 
in double-Weyl semimetals at low energy which implies 
that various measurable physical quantities receive loga¬ 
rithmic corrections at low temperature or low energy. 

Since the double-Weyl fermions disperse with different 
power behaviors along different momentum directions, 
the scaling dimensions of spatial coordinates can be dif¬ 
ferent. So, without loss of generality, we assume that 
scaling dimensions are given by [ 0 : 3 ] = —Z 3 ,[xi,X 2 ] = 
—zi,[t] = —1. In the following RG analysis, we inte¬ 
grate out the degree of freedoms in a momentum shell 
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FIG. 2. (a) Feynman diagrams, (b) The RG flow in the pa¬ 
rameter space of a = and j3 = ■ There are two 

fixed points. The fixed point at (0,0) represents noninteract¬ 
ing double-Weyl fermions which are not stable against long- 
range Goulomb interactions. The nontrivial fixed point is at 
(a*,/3*) = (0,1), where the Goulomb interactions at low en¬ 
ergy are qualitatively different between the quadratic- and 
linear-dispersing directions. 


oi h.e < A, where ^ > 0 is the flow pa¬ 

rameter, A is the ultraviolet energy cutoff, and = 

\j + v'^k'l which is the energy of double-Weyl 

fermions. We emphasize that momentum shell of this 
kind correctly includes the high-energy degrees of free¬ 
dom which effect the low-energy physics. From the var¬ 
ious coupling constants in the action, we construct two 
dimensionless coupling constants that will be useful to 
characterize the RG flow: a = and j3 = Note 

that the dimensionless parameter a, which has the form 
of a fine-structure constant, characterizes the strength of 
Coulomb interactions. 

We outline the main results of RG analysis here and 
provide calculational details in Appendix. We choose Zi 
and Z 3 such that the physical quantities m and v charac¬ 
terizing the dispersion of double-Weyl fermions are finite 
at the fixed point. Then, we can use only two parameters 
a and /3 to depict the full RG flow. By iterating the RG 
equations, we obtain the RG flow of a and j3, as shown in 
Fig. I^b). It is clear that there are only two fixed points: 
(a*,/?*) = (0,0) is the unstable fixed point of noninter¬ 
acting double-Weyl fermions while (a*, /3*) = (0,1) is the 
stable one. Near the fixed points, the RG equations to 
the lowest order of a and (3 are given by 


da 

m 

U 



(l-/3-ba)/3, 


( 4 ) 

( 5 ) 


from which it is apparent that the Coulomb interaction 
characterized by a is marginally irrelevant close to the 


stable fixed point of (0,1). It is known that marginally 
irrelevant interactions will give rise to logarithmic correc¬ 
tions to physical quantities in low energy or low tempera¬ 
ture. Besides being marginally irrelevant. Coulomb inter¬ 
action becomes infinitely anisotropic due to the screen¬ 
ing of anisotropic gapless double-Weyl fermions, namely, 
ryocQ://3—^-Oin the infrared. This phenomenon is cap¬ 
tured by the particle-hole polarization, 

16 

pi + VpI - n(p) ^pi{l + -^al) -b r]pl{l + pi). (6) 

At the stable fixed point (a*,/?*) = (0,1), we find that 
boson kinetic energy in the xy plane gets no qualitative 
renormalization, while it is strongly renormalized along 
the z axis. This renormalization can be understood as 
anomalous dimensions, 

f - n(p) +Py +P^[1 + ln^^] 

v\pz\ 

pI+pI + \pA- ( 7 ) 

In other words, the Coulomb potential is renormalized 
to be Vc{p) = p 2 _)_p 2 ^|p I ■ This result is qualitatively 
consistent with the one from the RPA calculation in Eq. 
([^, which we believe further shows that the divergent 
anisotropic screening of the Coulomb interaction is cor¬ 
rect. 

Screening of a Coulomb impurity: We expect that 
the anisotropic dispersion of double-Weyl fermions can 
cause the unusual feature of induced charge density of a 
Coulomb impurity, which interacts with the double-Weyl 
fermions through usual Coulomb potential proportional 
to 1/r. We shall consider both the case of noninter¬ 
acting double-Weyl fermions and of Coulomb-interacting 
double-Weyl fermions. For noninteracting double-Weyl 
fermions, a Coulomb impurity causes a divergent distri¬ 
bution of density of induced charges, while for Coulomb¬ 
interacting double-Weyl fermions, the induced charge dis¬ 
tribution is instead a simple 5 function located around 
the impurity. This 5 function feature of induced charge 
density is caused by effective screening of the Coulomb 
interactions by double-Weyl fermions. 

To be specific, we consider an impurity of charge Ze 

located at the origin, where Z is an integer constant. 

For noninteracting double-Weyl fermions, this impurity 

generates a bare Coulomb potential to them, namely, 
2 

Vo(^ = il^, where characterizes the Coulomb in- 
teraction between the impurity and the double-Weyl 
fermions. Then, the density of induced charges is given 
by Poij) = J where poiq) = ZeVoiq)Tl{u} = 

0 ,q), with Ii{LO,q) being the particle-hole polarization 
which we assume, without loss of generality, to be 
n(( 7 ) = — Bz\qz\ according to RPA calcula¬ 

tions. The anomalous characteristic of charge distribu¬ 
tion can be seen from partially induced charge densities, 
i.e., Qz{z) = f d?r_Lp{r) and Qxirx) = f dzpif) with 











4 


Q_L(r_L) cx (5(^)(r_L) and Qz{z) cx log|z| + const. Note 
that partially induced charge density Qz above is diver¬ 
gent at the long-wavelength limit, which arises because of 
the particular dispersion of noninteracting double-Weyl 
fermions. Such singularity does not appear when we con¬ 
sider Coulomb interactions in the double-Weyl fermions. 

For the interacting case, the Coulomb interaction is 
renormalized as V{^~^ = + q^ + |( 7 z|, where we have 

rescaled momentum to make the parameters look sim¬ 
pler. The induced charge density is now given by 

p{cf) = ZeV{^U{uj = 0,q} 

Ze{BA_qj_+ Bz\qz\) 
ql + q^y + \q.\ ’ 

from which we can compute Qz{z) and Q_L(r_L). Fur¬ 
thermore, / dzQziz) = J cPr±Q±{r±) requires that 
B± = Bz (see Appendix). Consequently, we obtain 
p{q) = BzZe, where Bz is a constant of order one. It is 
now clear that p{f) = BzZeS{r), where d{f) is the usual 
S function. This feature is similar to the induced charge 
discussed in Ref. [31], though the Coulomb interaction 
there is irrelevant. 

Logarithmic corrections to specific heat: Since 
the Coulomb interaction is marginally irrelevant at the 
nontrivial stable fixed point of ( 0 , 1 ), measurable physi¬ 
cal quantities will receive logarithmic corrections, which 
originate from the logarithmic behavior of the coupling 
constant under RG flow. More specifically, consider the 
fine-structure constant as a function of energy scale E, 


a{E) = 


oq 


1 + 


19q!q 

3 


log 


A > 
E 


(9) 


where Oq is the bare fine-structure constant and A is 
the ultraviolet energy cutoff. It is clear that as the en¬ 
ergy/temperature decrease, the interaction strength be¬ 
comes weaker in a logarithmic manner. One would ex¬ 
pect that this feature shall bring logarithmic corrections 
to various physical quantities. 

Near the fixed point, the free-energy density scales as 
/(T) = where 6 = ^ and zi, Z 3 are scal¬ 

ing dimensions of the xy coordinates and z coordinate, 
respectively [55] . For the specific heat C{T) = —T , we 
obtain C'(T) = Remember that we have 

set the scaling dimension of time to be one, while allow¬ 
ing zi,Z 3 to flow. As discussed earlier, close to the fixed 

point, we choose Zi = 5(1 — ciaydjS — ^J^ahia) and 
2 

Z 3 = 1 -b ^ where ci « 0.355, such that the mass 

m and Fermi velocity v in the dispersion of double-Weyl 
fermions remain fixed. Taking account of the logarithmic 
flow of fine-structure constant a, the specific heat at low 


temperature (for details see Appendix) is given by 

T2 

C{T) 


( 1 -^ i^ln A)5a 


X exp 


^ , ln(I + —In-) 
76V7r V 3 


, ( 10 ) 


where C 2 = ci — - « 0.037. The first line is the usual 

TT 

logarithmic correction to specific heat found in much of 


the literature. However, the second line in Eq. (10 1 , 


which actually dominates at low energy, is a different 
phenomenon that we find for double-Weyl fermions with 
marginally irrelevant Coulomb interactions. It decays 
faster than the usual logarithmic correction, which indi¬ 
cates that this marginal phenomenon deviates from Fermi 
liquid in a more radical way and may be tested in fu¬ 
ture measurements of HgCr 2 Se 4 and SrSi 2 at low enough 
temperature. It is crucial for the Fermi level to be at 
the double-Weyl points to observe the above logarithmic 
behavior of specific heat (see Appendix for a schematic 
specific heat curve). An interesting specific heat mea¬ 
surements in HgCr 2 Se 4 was reported recently [3^. How¬ 
ever, because the Fermi level in the measured material 
deviates from double-Weyl points, electron contribution 
to the specific heat does not show the logarithmic correc¬ 
tion in Eq. (10 1 . To possibly see the logarithmic correc¬ 


tions, it is desired to experimentally tune the Eermi level 
to the double-Weyl points in these candidate materials 
in the future. As discussed in Appendix, this unusual 
logarithmic correction is caused by the logarithmic term 
in the RG equations close to the fixed points. Logarith¬ 
mic terms also appear in RG equations describing cer¬ 
tain quantum criticality such as nematic quantum critical 
points in nodal d-wave suDerconductors [371 133] . 

Discussions: So far, we have considered a mini¬ 
mal model with only two double-Weyl points, e.g. in 
HgCr 2 Se 4 , interacting through Goulomb interactions. 
For such systems including HgGr 2 Se 4 , the Goulomb in¬ 
teractions are screened anisotropically: the screened 
Coulomb potential is qualitatively unchanged along the 
z axis, i.e., Vc( 0 , z) oc but behaves as Vc(’’-L) 0 ) oc ^ 
in the xy plane, which decays qualitatively faster because 
of screening by the low-lying quadratic fermions in the 
xy plane. In a material that possesses many double-Weyl 
nodes along various directions, e.g., in the predicted can¬ 
didate material SrSi 2 , one needs to take into account 
the contributions to particle-hole polarization from all 
double-Weyl fermions while the interactions between dif¬ 
ferent double-Weyl fermions can still be neglected as long 
as we care about low-energy and long-wavelength physics, 
namely, n(p) = ni(p), where i ranges over all double- 
Weyl points. In general, there is no special spatial direc¬ 
tion when many double-Weyl points are present. Thus, 
by averaging, we expect n(p) oc |^. In other words, the 
screened Coulomb potential has the form Vc(p) ~ 
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where |p| = +Py + p1 - In real space, it behaves like 

Vc(r) ~ ^ -vyhere Ji[x) is a Bessel function of the 

first kind and a is of the order of the lattice constant. 
For X » 1, Ji(a;) ^ Thus, the Coulomb potential is 

modified to be Vc{r) ^ 

To conclude, in this work, we explored the fate of 
double-Weyl semimetals in the presence of long-range 
Coulomb interaction. Instead of being screened to be 
short ranged, the Coulomb interaction is marginally ir¬ 
relevant at low energy. At the infrared stable fixed point, 
the Coulomb potential is strongly anisotropic, which 
is qualitatively different from the noninteracting Gaus¬ 
sian fixed point. Because of the marginal irrelevance of 
Coulomb interactions at the stable fixed point, physical 


quantities received logarithmic corrections in low energy 
or in low temperature. For instance, we find that the 
specific heat features an exotic logarithmic correction, 
deviating from the usual Fermi liquid, which originates 
from the interplay of anisotropic double-Weyl fermions 
and the bare isotropic interactions between them. 
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interesting work on a similar topic [35] ■ The RG scheme 
employed in our paper is different from Ref. [39] and some 
results are qualitatively different. 


APPENDIX 

1. Hubbard-Stratonovich transformation for Coulomb interaction 


Using Hubbard-Stratonovich transformation, we introduce a boson field to decouple Coulomb interaction. In double- 
Weyl semimetal, low-energy excitations are described by gapless fermions around the two double-Weyl points. Thus, 
with a momentum cutoff near band-touching points, we can write fermion operators as 


= 


d^k 

d^k 


4'(fc)e 




d^k 




ikr —iQr 


= -h V'-Q(a;)e 


-iQr 


(Al) 


where the electron’s orbital indices are implicit and 4>±Q(fc) = '^{k ± Q). Then, the charge density operator can be 
expressed by these low energy fermions: 


p{f) = 4'I(r’)4'(r) 
= PQi^+ 


(A2) 


where p±q{N) represent contribution from double-Weyl fermions around band-touching points. There are also mixing 
terms with nonzero momentum. Then, the (Euclidean) action for the instantaneous Coulomb interaction between 
electrons is given by 



d?rd?r' 


P{^P{^) 

47r|r — f'l 


d^rd^r' 


[pq(T) + P-Q{f)\[pQ{f) -f P-q(F')] 

47r|r — "PI 


-|- finite momentum terms. 


(A3) 


where = e^/e describes the strength of Coulomb interaction (e is the dielectric constant of the material in question). 
Note that finite momentum terms behave like f d^r ^ ; for a generic finite momentum Q, these terms are not 

important thanks to the fast oscillation factors which are neglected in our calculations. Using a Hubbard- 

Stratonovich transformation, we introduce a boson field ^(x) to decouple the Coulomb interaction, namely. 




Dtpexp 





(A4) 


We emphasize that boson field is not dynamic since it describes an instantaneous Coulomb interaction. This Coulomb 
potential plus a band Hamiltonian for double-Weyl fermion gives the effective action mentioned in main text. 
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2. Random phase approximation 


Coulomb interactions can be screened to some intent by the low-energy fermions near double-Weyl points. The 
screening effect by low-energy electrons is given by the boson’s self-energy appeared in RG analysis. Implementing 
random phase approximation (RPA), we shall consider only one loop diagram, which is first order quantum correction 
to boson self-energy (see next section): 

n(p)=5"/ ^Tr[roG(A:)roG(fc+p)] 

_a 2 f -w{w + rt) + {l/2m)'^{di{k)di{k+p) + d2{k)d2ik+p)) + v'^k3{k3+p3) 

^ Jk{w^ + i^/‘2m)^dlik)+dl{k))+v^kl){iw + il)^ + {l/2m)^dl{k + p) + dUk+p))+v^k3+P3)^)^ ’ 

where p = (w,p), J), = J and di{k) = kf — k"^, d 2 {k) = 2 kik 2 . At this rate, we employ Feynman parameters 
and dimensional regularization to evaluate this integral (rescaling the momentum to absorb coupling constants) 


22 m f -w{w + ft) + k 3 {k 3 + P 3 ) + di{k)di{k + p) + d 2 {k)d 2 {k + p) 

[p)- 9 — + {kl + klY)[{w + ny + {k3 + P3? + ((fci + piY + (fc 2 + P 2 m 

_ 2 g^ 2 m f —x{l — x){—i}'^+pl)+di{k)di{k+p) + d 2 {k)d 2 (k+p) 

27r V Jx,k± x(l — x){fl‘^ + pI) + {1 — x)kj_ + x\k± + pxl'^ 


where x is the Feynman parameter, dx, and k± = (^ 1 ,^ 2 )- It is clear that n(n = 0,p = 0) = 0 which is 

correct because of the vanishing density of states at the double-Weyl points. In the limit of low-energy, the momentum 
transfer is small so that we can obtain the polarization approximately. Momentum dependence of self-energy along 
P 3 direction are screened by the quadratic dispersion of fermion spectrum: 


n(fl = 0 ,pj_ = 0 ,p 3 ) = 


The dependence on px and Py is given by: 

n(fl = 0,pj_,p3 = 0) = ^ — 

ZTT V . 


(X -pj_, 


2 g^ 2 m 
2 tt V 


—x{l — x)pI + k\ 

X{l-x)pl + kj_ 


2 g^ 2 m 
2 'k V 

oc -|P3|- 


-1 


P3 


8 k'\ tanh , - 

-L V4fel-epi 


P3 




-pI 


(A7) 


’ x,kj 


{kj - fc|)((fci -I-Pl)^ - (fe +P 2 )^) -|-4fcifc2(fcl +Pl)(fc 2 +P 2 ) 
{1 - x)kj_ + x{k± + p±)‘^ 


(A8) 


where we have implicitly neglected a potentially-unimportant logarithmic factor logpj_. Consequently, we obtain 
electron-hole polarization: 

n{p) (X-{p‘i + \p3\). (A9) 


We note that the momentum dependence of Coulomb interaction along the P 3 direction is strongly screened by electron- 
hole polarizations. Since at low energy, linear behavior of momentum p 3 is more important than bare quadratic 
momentum dependence, this will dramatically affect the Coulomb interaction: it becomes infinitely anisotropic. This 
anisotropy can be seen more directly in real space. The Coulomb interaction along z axis is like usual R(0,a;3) oc |^; 
however, in the xy plane, it behaves as V{0,x,y) oc yrpp-. The effective interaction in the xy plane after screening is 
qualitatively weaker because quadratic dispersions in the xj/-directions render larger density of states at low energy 
and can screen the interactions more effectively. 


3. Renormalization group analysis 


Here we perform renormalization group (RG) analysis in details, 
double-Weyl fermions with long-range Coulomb interactions is 


S= d'^x < Ip 


{do+ig(l))ro - —dj{-iW)rj - iv{-id3)r3 
2 m 


The effective (Euclidean) action describing two 
^ + {d 2 <pf) + v{d3(l)T)'^ , (AlO) 
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where the ultraviolet cutoff is implicit in the above equation. We can construct two dimensionless parameters from 

2 2 

various coupling constants in the action above a = and /3 = . In Wilsonian RG, integrating out the high- 

energy modes with /r < E’g < A, where will generate an effective action with lower 

energy cutoff. We derive the RG flow from iteratively integrating momentum shell whose fermions having energy 
G (Ae“*, A), where ^ > 0 is the RG running parameter. At the one-loop level, the effective action reads 

= + (All) 

where Sc is the boson-fermion coupling which involves high energy modes in the momentum shell and (• • •) means 
taking expectation value from high energy part of field configurations. There are three relevant contributions, i.e., 
boson self-energy, fermion self-energy and vertex correction. It is worth noting that since the boson field is not 
dynamical, fermions do not get an anomalous dimension. Moreover, vertex correction at one-loop level vanishes. This 
is related to Ward identity. 

It is straightforward to obtain particle-hole polarization and fermion self-energy through direct calculations as 
follows: 


n(fc) =g^J ^TT[roG{q)ToG{q + k)] 

2g^l ki+k2 g^mvl rj 2 

2 SttAt] 2 

nk) = -g^ j ^T^G{q)T^D{k-q) 

g‘^xFi{x)l , ^ g^xE2ix)l I ^,,2 , 2 ^T. , 07 7 t-i 

— -vk^Ts --— [(fci - fc2)ri + 2fcifc2r2j , 


Snv 


8ttv 2m 


(A12) 


(A13) 


where x = and Ei,E 2 are two positive functions whose behaviors will be discussed below. Including boson and 
fermion self-energy, we obtain the full effective action at lower energy cutoff: 


SeS = S<+ + i ^ dV<(-n)<(><. 

After rescaling, we could obtain RG equations of all coupling constants: 

=2zl-l-G2{x)^/^, 


dl 

dhiv 
dl 

d\ng^ 8 1 

— = l-«-3a--A 


= 1- Z 3 + Gi{x)a, 
16 


(A14) 

(A15) 

(Alb) 

(A17) 

(A18) 


where x = and Gi{x) = x'^Fi{x‘^),G 2 {x) = x^F 2 {x‘^) are two well-defined positive functions which 

are plotted in Fig. Here, Zi is the scaling dimensions of spatial coordinates Xi which were defined in main 
text. We assume that the fermion mass and Fermi velocity are finite at the fixed point since they are measurable 
quantities defining the dispersion of fermions. Consequently, the scaling dimensions of spatial coordinates should 
satisfy Z 3 = 1 -I- Gia and Zi = 5(1 + G 2 Vctp). Note that fermions are free from anomalous dimension due to 
instantaneous Coulomb interaction. The RG flow is best captured by the RG equations of dimensionless coupling 
constants a = and /3 = f^: 


dcy, 2r / \ ^^1 

- = -a |,G.(x) + -], 

d^ 


f5[l- P + a-xGi^x) - v^7ra/3G2(a;)], 


(A19) 

(A20) 
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Gi(x) G 2 (x) 



FIG. 1. The numerical values of the Gi functions are plotted. 


where the functions Gi and G 2 are well-defined on real axis. Moreover, we know their values at zero and infinity: 
Gi(0) = 02 ( 0 ) = 0 and G'i(oo) = 0, G' 2 (oo) = Near x = 0, they can be Taylor expanded as 


Gi(x) = --^ + 0(x3), (A21) 

TT 4 

G 2 {x) = —(—3 -I- 41n2 — 41nx)x -I- 0{x^). (A22) 


By iterating above coupled RG equations (see the RG flow diagram in the main text), it turns out that there are 
two fixed points with {a*, P*) = (0,0) and (a*,/3*) = (0,1). The former represents the noninteracting fixed point 
and is unstable. However, the latter is a stable fixed point. It is worth to point out that r] flows to zero at infrared 
even though the bare value of 779 is close to one because of isotropic Goulomb interaction in the ultraviolet. This is 
qualitatively consistent with RPA calculations that low-energy double-Weyl fermions screen the Coulomb interactions 
in an anisotropic way. For simplicity, we employ the following linearized RG equations to capture the two fixed points: 


da 19 2 

^ = (l-/3 + a)/3, 


(A23) 

(A24) 


from which it is clear that (a*,/3*) = (0,1) is stable fixed point while the noninteracting one (a*,/3*) = (0,0) is 

unstable. At the fixed point of noninteracting double-Weyl fermons, /3 is strongly relevant while a is marginally 

irrelevant. Relevant /3 around the fixed point of (0,0) drives the system away from this unstable fixed point. On the 

2 

other hand, near the stable fixed point (0,1), a is marginally irrelevant while /3 is irrelevant. Since a oc ^ describes 
the strength of Coulomb interaction, we conclude that Coulomb interaction is marginally irrelevant in the infrared 
limit, which gives rise to logarithmic corrections to physical quantities as we mentioned in the main text. For the 
stable hxed point, because a* = 0 the scaling properties of spatial coordinates are same at the tree level, namely, 
Zi = ^, Z 3 = 1 at the fixed point. 


4. Screening of a charge impurity 

Consider the case of noninteracting double-Weyl fermions first. As mentioned in main text, induced charge density 
in momentum space by a Coulomb impurity is given by 


p(g) = ZeR(g)n(c. = 0,<7) 

2 B±q\ + B3I93I 


— 


dl + ql 


(A25) 












9 


from which we can compute the density of induced charge in the real space by the Fourier transformation: p{f^ = 
I denote p(r) = pi(fj + pii(r^, where 

1 


Pi(r±,r3) = -ZegoB± 


= -ZeglB± 


= -ZeglBj_ 
= -ZeglBj_ 


= -Zeg^B^ 


{2lTf 

TT 

TT 

TT 

(27r)2 


J r± + % 

/ dg_Lgie"l’'3l9i 


2rQ 


(27r)2 (r^ + r|)5/2 ’ 


and 


1 


d q±dq3 


ksl 


<71 + 


7 2 0 /■ .2 ^2,1 f q\rl 

= -ZegoB3j^ j d — 


o-^Q3^3-'^Q±'>^± 


0 0-*^ 
^7 '-'7 2 

2 ^2 


-tq±r± 


^®^°^^(27r)2 

Put them together, we have induced charge density: 

p(f) = pj(f) +p„(r) 

1 / „ —r 


-^^9oB3-j0^ J dq±q±Jo{q±\r\)Glll 

1 2(^r2_ + - r3 sinh"^(T^)) 


0,0, i 


(r^ + r|)3/2 


i+2ri 


= (2^ r^^ (ri+r|)^/2 + 


yriTrf-rssinh 


(r^ + r|)3/2 


Partially induced charge densities introduced in the main text are given by 

Q±ir±) = J dr3p{rj_,r3) 


= -Zegl J dr3 J 

i2nyql' 

= -ZeglB^S^^\r^), 


dq3 dql BjLql + B3\q3\ 


2tt (27r)2 qj_ + g| 

dql ql-i„^.r^ 


Qsirs) = J d^rj_p{rj_,r3) 
= -Zegl J d^rj_ J 


d^ dq\ B3^q\ + B 3 \q 3 \ 


2tt (27r)2 ql_ + gf 


= -ZeglB3 


dq3 1 
27r 1931' 


^q3r3 




(A26) 


(A27) 


(A28) 


(A29) 


(A30) 


where 7 is Euler-Gamma constant. 

We now consider the case of interacting double-Weyl fermions. In a similar way, the induced charge density in 
momentum space is given by 

Pid) = Zeg‘^V{q)n{uj = 0,q) 

2B±ql + S 3 1931 


= -Zeg^ 


ql + I<73| 


(A31) 
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from which we obtain partially-integrated induced charge density in real space: 


Q±{r±) 

Qsirs) 


- / drs 


^ ^ (27r)3 

-Zeg^B^S^^\r^), 

-]d\^ I 

-Zeg^BsS^rs). 


^^^2 B±ql + Bslqsl 


dl 


1931 


(2^)3 ^ ql + |g3| 


(A32) 


(A33) 


As the total induced charge 
is required that B± = B^. 
-Zeg^B3S^^\r^. 


should not depend on the order of integration, i.e., J d^r±Q±{r±) = f dr 3 Q 3 {r 3 ), it 
So, induced charge density in real space is p{f) = —Zeg'^B^ f = 


5. Specific heat 


Since Coulomb interaction is marginally irrelevant at the stable fixed point, we expect logarithmic corrections to 
physical quantities. Here, we specifically calculate how specific heat of double-Weyl fermions are affected by Coulomb 
interactions in low temperature. Corrections to other measurable quantities may be obtained similarly by taking 
account of scaling hypothesis near the fixed point as well as the flow of coupling constants. According to scaling 
hypothesis, free energy density scales as / = where b = ^ and zi,Z 3 are scaling exponent of coordinates 

along X (or y) and 0 directions respectively. As a result, specific heat C = —T is given by 

CiT) = &-2^i-^^C'o(A), (A34) 


from which we obtain the RG equation of specific heat as follows: 

^ = -(2zi + Z3)C'(0, (A35) 

di 

where I = \nb. Near the fixed point, zi = 1(1-1- G 2 \/^) and Z 3 = 1 -I- Gia (see the preceding section of this 
Appendix for details). Because a is magically irrelevant and /3 is irrelevant near the stable fixed point, the RG flow 
along P direction is much faster than that along a direction. This is obviously seen in Fig. 2(b) of the main text. As 
logarithmic correction is caused by marginal flow of a, we can set /3 = /3* = 1 in the following calculations as a good 
approximation. For a close to a* = 0, we obtain the asymptotic behavior of Gi and G 2 : 


Gi 


1 

TT 


TT 



G 2 ~ —ci\/a — 



(A36) 

(A37) 


where Ci > ^ is a positive constant. Note that the logarithmic term in G 2 is due to fermion self-energy. By plugging 
Eq. (A36), Eq. (A37l, and Eq. (9) of the main text into Eq. (A35), we obtain: 


dlnG 

dl 


— —2 — C2<y.{pj -\- 


•y/Hr 


a{l) Ina(^), 


(A38) 


where a{l) = —^ 5 - and C 2 = ci — - « 0.037. The first term will bring the usual logarithmic correction to specific 

heat as we expect. However, the second term becomes more dominant in the infrared or at low enough temperature. 
The deviation of specific heat from the usual Fermi liquid behavior is apparent: 


G{T) oc exp 


19 ^ 


T2 


19ao 1 ^0 3 


(1 + 


X exp 


760? 


ln(l 


ln(l + 0-1^0) 


19ao , Tq ^ 

—‘"r> 


(A39) 


The exponential factor in the last line of Eq. (A39) is a new phenomenon of marginal Fermi liquids we find in double- 
Weyl fermion (the usual logarithmic corrections appear without the second exponential factor in Eq. (A39)). Fig. 
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c/r^ 



FIG. 2. Schematic plot of specific heat in double-Weyl semimetals with and without Coulomb interactions. The dotted l ine is 
the specific heat from electronic contribution in noninteracting double-Weyl fermions. The solid line is plotted from Eq. (A39l 
[or Eq. (10) in the main text]. The dashed line is plotted from the same equation but without the exponential factor. 


(S2) shows the schematic plot of specific heat. The solid line is plotted from Eq. (A39l and is more exotic comparing 
to the dashed line which is the usual logarithmic corrections appearing near common marginally-irrelevant fixed point. 
The origin of this more exotic corrections lies in the unusual logarithmic contribution in G 2 , which comes from the 
interplay between the anisotropic spectrum of double-Weyl fermions and isotropic bare Coulomb interactions. 


* yaohong@tsinghua.edu.cn 

[1] O. Vafek and A. Vishwanath, Annu. Rev. Condens. Mat¬ 
ter Phys. 5, 83 (2014); and references therein. 

[2] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim Rev. Mod. Phys. 81, 109 
(2009). 

[3] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ). 

[4] X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 

( 2011 ). 

[5] S.-S. Lee, Phys. Rev. B 76, 075103 (2007). 

[6] Y. Yu and K. Yang, Phys. Rev. Lett. 100, 090404 (2008). 

[7] T. Grover, D. N. Sheng, and A. Vishwanath, Science 344, 
280 (2014). 

[8] S.-K. Jian, Y.-F. Jiang, and H. Yao, Phys. Rev. Lett. 
114, 237001 (2015). 

[9] K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. 
Rev. Lett. 103, 046811 (2009). 

[10] Q. Liu, H. Yao, and T. Ma, Phys. Rev. B, 82, 045102 

( 2010 ). 

[11] O. Vafek and K. Yang, Phys. Rev. B 81, 041401(R) 

( 2010 ). 

[12] Y. Lemonik, 1. L. Aleiner, C. Toke, and V. 1. Falko, Phys. 
Rev. B 82, 201408(R) (2010). 

[13] W.-F. Tsai, C. Fang, H. Yao, and J. Hu, New J. Phys. 
17, 055016 (2015). 

[14] A. A. Abrikosov, Sov. Phys. JETP 39, 709 (1974). 

[15] E-G. Moon, C. Xu, Y. B. Kim, L. Balents, Phys. Rev. 
Lett. Ill, 206401 (2013). 

[16] 1. P. Herbut and L. Janssen, Phys. Rev. Lett. 113, 106401 
(2014). 

[17] H.-H. Lai, B. Roy, and P. Goswami, arXiv: 1409.8675 

[18] P. Goswami and S. Chakravarty, Phys. Rev. Lett. 107, 
196803 (2011). 


[19] P. Hosur, S. A. Parameswaran, and A. Vishwanath, Phys. 
Rev. Lett. 108, 046602 (2012). 

[20] H. Isobe and N. Nagaosa, Phys. Rev. B 86, 165127 

( 2012 ). 

[21] H. Isobe and N. Nagaosa, Phys. Rev. B 87, 205138 (2013). 

[22] W. Witczak-Krempa, M. Knap, and D. Abanin, Phys. 
Rev. Lett. 113, 136402 (2014). 

[23] H.-Z. Lu, and S.-Q. Shen, Phys. Rev. B 92, 035203 
(2015). 

[24] X. Wan, A. M. Turner, A. Vishwanath, and S.Y. 
Savrasov, Phys. Rev. B 83, 205101 (2011). 

[25] A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 
127205 (2011). 

[26] W. Witczak-Krempa and Y. B. Kim, Phys. Rev. B 85, 
045124 (2012). 

[27] H. Weng, C. Fang, Z. Fang, B. A. Bernevig, and X. Dai, 
Phys. Rev. X 5, 011029 (2015). 

[28] B. Q. Lv, H. M. Weng, B. B. Fu, X. P. Wang, H. Miao, J. 
Ma, P. Richard, X. C. Huang, L. X. Zhao, G. F. Chen, Z. 
Fang, X. Dai, T. Qian, and H. Ding, arXiv: 1502.04684 

[29] S. Y. Xu, I. Belopolski, N. Alidoust, M. Neupane, C. 
Zhang, R. Sankar, S.M. Huang, C.C. Lee, G. Chang, B.K. 
Wang, G. Bian, H. Zheng, D. S. Sanchez, F. Chou, H. 
Lin, S. Jia, and M.Z. Hasan, arXiv: 1502.03807 

]30] G. Xu, H. Weng, Z. Wang, X. Dai, and Z. Fang, Phys. 
Rev. Lett. 107, 186806 (2011). 

[31] S.-M. Huang, S.-Y. Xu, 1. Belopolski, C.-C. Lee, G. 
Chang, B. K. Wang, N. Alidoust, M. Neupane, H. Zheng, 
D. Sanchez, A. Bansil, G. Bian, H. Lin, and M. Zahid 
Hasan, arXiv:1503.05868 

[32] C. Fang, M. J. Gilbert, X. Dai, and B. A. Bernevig, Phys. 
Rev. Lett. 108, 266802 (2012). 

[33] H. B. Nielsen and M. Ninomiya, Phys. Lett. B 130, 389 
(1983). 

[34] B.-J. Yang, E.-G. Moon, H. Isobe and N. Nagaosa, Nat. 
Phys. 10, 774 (2014). 

[35] D. E. Sheehy and J. Schmalian, Phys. Rev. Lett. 99, 








12 


226803 (2007). 

[36] W. D. Wang, A. Li, T. Dong, M. Lei, X. L. Fu, S. S. 
Miao, P. Zheng, P. Wang, Y. G. Shi, J. L. Luo, and N. 
L. Wang, J. Low Temp. Phys. 171, 127 (2013). 


[37] Y. Huh and S. Sachdev, Phys. Rev. B 78, 064512 (2008) 

[38] E.-A. Kim, M. J. Lawler, P. Oreto, S. Sachdev, E. Frad- 
kin, and S. A. Kivelson, Phys. Rev. B 77, 184514 (2008) 

[39] H.-H. Lai, Phys. Rev. B 91, 235131 (2015). 



